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Abstract 

A toy model of strongly correlated fermions is studied using Green 
function and functional integration methods. The model exhibits a 
metal-insulator transition as the interaction is varied. In the case of 
unrestricted hopping is established an equivalence of the model with 
the Hubbard model with infinite range hopping. The generalization 
to the case with N components is made. 

1 Introduction 

The simplest model to describe strong electron correlations is the Hubbard 
model Unfortunately, the only exact solution available corresponds to the 
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one dimensional case . Thus, to understand its physical properties we have 
to relay upon approximate solutions which in many instances are highly non- 
trivial. Another alternative is to formulate exactly soluble toy models which 
share some properties with realistic many body models. The importance of 
the toy models lies in the fact that many interesting properties of complicated 
strongly correlated electron systems are simulated in a compreensive way. 

A simple toy model was proposed by Hatsugai and Kohmoto (to be 
referred from now on as the 'HK model'). This model is exactly solvable 
in a very simple way and have the same atomic and band limits as the 
Hubbard model. Also, as Hatsugai and Kohmotto have shown, the model 
describes a metal-insulator transition (MIT). This MIT scenario of the HK 
model was studied from a scaling point of view by Continentino |^ who also 
formulated a boson version of the HK model. A closely related toy model 
was also discussed by Baskaran |Q. There is also an atractive version of the 
HK model which was proposed recently by Mattis and Bedersky 0. 

Another different toy model consisting by a Hubbard model with an in- 
finte range hopping has been solved exactly ||^. This model has the same 
atomic limit as the Hubbard model but does not have obviously the same 
band limit since hopping is unrestricted. As shown in ref. [^, this model is 
an insulator at half-filling for any on site Coulombian repulsion U > 0. 

In this work we discuss the HK model using Green function and integral 
functional methods. It is shown that the HK model and the Hubbard model 
with infinite range hopping are equivalent in the thermodynamic limit if the 
hopping in the HK model is assumed to be unrestricted. We also solve exactly 
the component HK model for any finite N. This solution is obtained by 
an exact evaluation of the functional integral representation of the model. 
From this solution one obtains also the exact solution of the N component 
Hubbard model with infinite range hopping. For the sake of completeness, 
we also dicuss a lattice fermion model where the interaction, rather than the 
hopping is of infinite range. 

2 The HK model 

The HK model is described by the Hamiltonian: 
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H = EEt^^) - /^]%. + (2-1) 

fc k 

where nj:^=c^j^^Cf:^ and e{k) = — 2t X^^Li cos/cj. yU is the chemical potencial. 

The dispersion e{k) corresponds to a nearest neighbour hopping on a hyper- 
cubic lattice in d dimensions. Note that, in contrast to the Hubbard model, 
the Coulombian repulsion is local in fc-space rather than in real space. Since 
the Hamiltonian is local in /c-space and the interacting part commutes with 
the non-interacting one, it is straightforward to solve it. For instance, the 
exact free energy density is 

/ = -i- ^ln{l + 2e^[^-^(^^W + e^[2M-2e(A^)-c/]| (3.2) 
k 

where (3 = 1 /T, T being the temperature, and L is the number of lattice sites. 
Note that in the zero bandwidth limit we obtain the free energy density of 
the atomic limit of the Hubbard model. Also, in the limit f/ = we get the 
band limit. 

The Matsubara Green function is defined by G^{k, r) = — < rc^^(r)ct^(0) >. 
In frequency representation it is given straightforwardly by 

-» 1 — 11 ^ 11 ^ 
G.ik, u^n) = ^ + : ^^TTT^' (2-3) 

where Un = (2n + l)7rT, nEZ, and 

= I Zl (2 4) 

Note that we are assuming a paramagnetic phase, that is, n^^ = n^^^. At 
T = Eq.( |2.4| ) assumes the form of a step function: 



n,^^(T = 0) = ^(^0 - e{k))[-e{U - \e{k) - /io|) + e{\e{k) - /io| - f/)], (2.5) 

where /io is the chemical potential at T = 0. 9{x) is the usual Heaviside 
function. This occupation number has two discontinuous jumps. An imme- 
diate consequence of this is that Luttinger's theorem does not hold. Note 
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that the exact Green function for the HK model share the same atomic and 
band hmits with the Hubbard modeL 

As in the case of the Hubbard model, a half-filling condition is obtained 
by setting /i = U /2. At half-filling and T = we have that the Green 
function becomes 

G.{k.u:) = - 7^ (2.6) 

In above, lo is not more a Matsubara frequency but an Euclidian frequency, 
that is, we still have an analytic continuation to imaginary time but in zero 
temperature. Thus, the frequency uj is continuous and not discrete. 

From the poles of Eq. ( |2.6| ) we obtain two bands in straight analogy with 
the case of the Hubbard model where we have the so called lower Hubbard 
band and the upper Hubbard band. The poles gives the energy bands: 



E. = e{k)-- (2.7) 

E+ = e{k) + ^ (2.8) 

Note that the gap has a size which depends on U. If the two bands are 
separated by a nonzero gap we have an insulator. The value of U which 
gives a zero gap is determined by demanding that the top of the lower band 
coincides with the bottom of the lower one. Thus we find that the critical 
value of U that signals a MIT is given by Uc = 4:td = W. Thus for U < Uc 
the system is metallic. Note that in the HK model we have a MIT at any 
dimension while in the case of the Hubbard model it is well known that in 
one dimension no MIT exists at half-filling 0. 

It is also interesting to study the behavior of the spin suscetibility. This 
is given by the response function: 

X{k, Vn) = -/x|n(fc, !/„) (2.9) 

where /is is the Bohr magneton and Il{k, z/„) is the polarization Green func- 
tion which is given by: 

n(fc, i^n) = 7^ E E G{k + q, Vn + i^n)G{q, CJ„), (2.10) 
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where Un = 2n7rT, n^Z, is a Bose Matsubara frequency. We are omitting 
spin indices for simplicity. The Matsubara sum is over a fermion frequency. 
We will use the Green functions at half-filling. We evaluate the suscetibility 
at zero temperature in the static limit and at the nesting wave vector Q = 
(vr, ...,7r). Performing the Matsubara sum and taking the zero temperature 
limit gives 



m - -/i f ^ f eiU/2 + e{q))-eiU/2-eiq)) 
^ 4 JiBz(27rY 2e(a) 



+ 



4 JiBZ {2ttY \ 2t{q) 

e{U/2 + t{q))~e{-e{q)-U/2) e{t{q) - U /2) - OjU /2 - t{q)) 
2e{q) + U 2t{q) - U 

^ e{e{q)-U/2)-e{-e{q}-U/2) \ 



When U = Uc = W the integral of the second term in Eq.( p.ll|) is divergent 



for q = ("infrared divergence") while the integral of the third term is 
divergent for g = Q ("ultraviolet divergence"). 

Let us evaluate explicitly the spin susceptibility. If we use a square density 
of states, p{uj) = {l/W)e{W/2 - uj)9{uj + W/2), and consider U >UcWe get 

^(g,0) = -i^lnf^^) (2.12) 

Thus we have in fact that the susceptibility diverges as U approaches Uc- 
Note that we do not have a simple power law behavior for the suscetibility. 
This should be contrasted with other descriptions of the MIT. For exemple, 
the Brinkman and Rice transition gives a power law behavior for the 
static suscetibility, which is calculated for a zero value of the reciprocal lattice 
vector rather than in the nesting wave vector. 

The analysis of this section allows us to conclude that the MIT in the 
HK model is a second order transition with mass gap A = U — Uc- Note that 
A plays a role similar to the ro coupling associated to the quadratic part in 
the Ginsburg-Landau Hamiltonian. Remember that in a Ginsburg-Landau 
theory it is assumed that ro~T — where is the critical temperature 
associated to the second order transition. Here we have U rather than T. 
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3 The Hubbard model with infinite range hop- 
ping 

The kinetic energy of the HK model is identical to the one of the Hubbard 
model. Thus it is natural to ask about the differences of behavior between 
the HK model and the Hubbard model as we manipulate or change the 
hopping matrix elements. In this section we will consider the HK model 
with infinite range hopping (IRH). The IRH Hubbard model was already 
studied and solved exactly [0]. We will solve exactly the IRH HK model and 
show that the solution is the same in the thermodynamic limit, to that of the 
IRH Hubbard model. This establishes the equivalence between the IRH HK 
model and the IRH Hubbard model. In the case of infinite range hopping it 
does not matter if the electrons interact locally in /c-space or in real space. 
In both situations the solution is the same. 

When the hopping is of infinite range the dispersion is given by e{k) = 
—Lt6j:Q 0. The exact free energy density for the IRH HK model is obtained 
by substituting this dispersion in the exact expression Eq.( p.2|) and taking 
the thermodynamic limit. It is readily obtained that 

fiRH = - ^ + 2e^^ + 6^(2^^-^)], (3.1) 

which is the same expression that corresponds to the free energy density of 
the IRH Hubbard model 0. 

Also, by substituting the IRH dispersion in Eq.( |2.3D and taking the ther- 
modynamic limit one obtains the exact Green function for the IRH HK 
model: 

GlRHik, Un) = -. ^ — + ^ — , (3.2) 

Zq being the atomic partition function per site. Since the free energy of the 
IRH HK model coincides with that of the IRH Hubbard model, the above 
expression gives the exact Green function of the IRH Hubbard model. It 
turns out that the two models have the same perturbative expansion in the 
thermodynamic limit. Note that the k = mode does not propagate while all 
the others k modes have k independent propagators. Note also that the k^O 
modes have their propagators identical to the atomic limit Green function. 
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This behavior is manifested also in the expression for the free energy, Eq.( p.l| ), 
where we have a term given by the atomic solution corresponding to the 
modes and a term —2t corresponding to the = mode. This means 
that in the IRH regime the = and modes separate. 
From Eq. (|3.2| ) one obtains the electron density: 



n 



1 + 26^/^ + e^(2/.-c/) 



(3.3) 



The above relation can also be obtained directly from the expression for the 
free energy as = —{dfiRH/dfi)/3. Solving Eq. (|3.3|) for fi one obtains 



fi = U + p-Hii 



(1 - ny + n{2 - n)e-P^ - (1 - 



2-n 



(3.4) 



in agreement with reference |^. 

From Eq. (|3.2| ) we obtain the density of doubly occupied sites d: 



d{P,U,fx) 



,/3(2/.-C/) 



(3.5) 



1 + 2e/3^ + e/5(2/.-c/) 

The above expression in function of n is far more compact than that one 
obtained in ref. which has been written as a function of n. When fi = U/2, 
which corresponds to n = 1, one has 



d 



2(ef^ + l) 



(3.6) 



which again agrees with ref.0. 



4 The HK model with N components 

In section 2 we solved exactly the HK model with = 2s + l = 2 components 
(s is the spin). In this section we show that this model is also exactly soluble 
for an arbitrary number of components. 
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Let us define the component spinors: 



\ ^kN J 



(4.1) 



(4.2) 



Each component of the expression above is a Grassmann field. We will 
evaluate the exact partition function by writing it as an functional integral 
over the fields defined above. The partition function is given by 



(4.3) 



where the action 5" is 



k 

where it is understood that the fields are time dependent (Matsubara time) 
and that they satisfy antiperiodic boundary conditions in Matsubara time. 

We eliminate the quartic term by means of a Hubbard-Stratonovich trans- 
formation, getting the new action, 

S' = dr - - C//2 + e(^) - iU4>^)^^, + (C//2)0|] (4.5) 

k 

where (f)^ is an auxiliary Bose field. Since the new action is quadratic in the 
Fermi fields, it is straightforward integrate out these and obtain the following 
effective action: 



Seff = -N 



lndet[9^ - // - C//2 + e(^) - iUct)^] + - dTj2(l>l (4-6) 

2 t/ -* 
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We can evaluate exactly the determinant appearing in the above equation by 
solving the differential equation: 

U/2 + e{k) - tU<P^{T)]Uk, r) = a^Uk, r) (4.7) 
with fn satisfying the antiperiodic boundary condition fn{k, 0) = —fn{k,(3). 



The Eq. ([4.7|) is easy to solve and the solution is given by 



fnik, r) = cexp dr'lfx + U/2 - e{k) + 2f/0,-(r') + , (4.8) 

c being an arbitrary constant. By applying the antiperiodic boundary con- 
dition to this solution we get 

anik:) = -iun - /i - f//2 + e(fc) - tiU//3) f drcj)^ (4.9) 



where ujn is a Fermi Matsubara frequency. Thus, the determinant which 
appears in the effective action is given by the product of the a„'s. Therefore, 



+00 ^ rl3 U 1-0 

Seff = -N Yl Y.H-^^n-^^-U/2+e{k)-^{U/P) dT<p^,{T)] + - drct>l{r 

n=-oo j: JO Z JO 

(4.10) 

The Matsubara sum in Eq.( [4.10D is easily done by standard methods. The 



partition function is then written in the form: 



Z = J +e^[^+^/'-^('^)le*^/o''^-<^eW|'^^^ (4.11) 

By expanding the term between braces and performing the Gaussian inte- 
grals one obtains the exact expression for the partition function of the N 
component HK model: 

^ = n E „ e-^t^-(^^)ie^"^ (4.12) 

From Eq. (gl2D one obtains the free energy density: 
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/ = -— Vln ^ V ^, e"^'^-(^^)]e^^^ . (4.13) 



Note that in the case of = 2 the result given in Eq.(2.2) is recovered. Also, 
in the zero bandwidth limit one obtains the atomic limit of a Hubbard model 
with A^ components. 

It is worth to emphasize that the above result is exact for any finite A^. 
We can also look for the solution at large A^. The limit A^ ^ oo of Eq.( 4.13| ) 



is not obvious because of the factorials. However, we can obtain the large 
N solution by a saddle point evaluation of the functional integral with the 
effective action ( |4.1CI|) . In order to perform this saddle point evaluation it 



is useful to make the rescaling U ^ j^, (p^^Ncj)^. In this way we get that 
5*6// = NSeff where Seff is A^ independent (note that the term U/2 inside 
the logarithm is rescaled to U/{2N) which is zero at large A^). In this case 
the saddle point solution corresponds to the exact solution when N ^ oo. 
The large A^ solution is found to be just the Hartree-Fock approximation for 
the HK model. 

From the exact expression ( 4.13| ) it is possible to obtain the exact free 



energy density for the A^ component IRH Hubbard model. It is given by 



fiRH{N) = -Nt-^\n 



N 



ton\{N-n)\ 



(4.14) 



5 The infinite range interaction case 



Let us write the interaction of the HK model in the real space representation. 
This is achieved by Fourier transformation: 



—7=J2^wi-ik-Ri)cia. (5.1) 

It is readily obtained that the interaction part is given by 

U 

ijkl 



Hi = —Y. ^i+k,j+ic^Cj^cl^cii. (5.2) 
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The above interaction is obviously of infinite range. We note that it is not 
necessary to consider only interaction between opposite spins. Generalizing 
the above interaction in a such a way as to include interaction between like 
spins correspond to a shift in the chemical potential. In this situation we have 
a particular case of the HK model which consist of a model with Hamiltonian 
given by 

<i,j> T i a i,j 

where rii — Z^f^rijo-- Since the interaction is the same for every pair of sites 
in the lattice we can write ntrij — (X^i ni)'^. It is straightforward to write 

— * 

the above Hamiltonian in /c-space as 

H = Ei:m - l^]n^,^ + ^(E^)'- (5-4) 

k k 

We write the partition function in the same way as in the last section, 
that is, by going to an integral functional representation, the fermionic fields 
being substituted by Grassmann fields. As before, the quartic term in the 
action is eliminated by applying a Hubbard- Stratonovich transformation. We 
obtain the following action involving an auxiliary Bose field 0: 



[ drii: E i^kaidr - + e{k) - iU<t>]i^^,^ + —0^}. (5.5) 
JO ^ „ I 

By integrating the fcrmions and carrying the calculations in a way similar to 
the one in the preceeding section, we get an effective action S'e// = LSeff, 
where Seff is given by 



— u fP U fP 

Seff^-2 depo{e)^]n[-iuJn- fJ^ + e-i— drcpir)] + — dT(f)'^{T), 

J—oo p Jo Z Jo 

(5.6) 

where po(f) is the bare density of states. Therefore, we write the partition 
function in the following form: 

Z ^ N [ D(pe-^^^^f^'^\ (5.7) 
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where is a normalization factor. From Eq.( |5.7| ) we see that in the ther- 
modynamic hmit the saddle point of the action corresponds to the exact 
solution of the problem. The saddle point is given by an imaginary number 
in the repulsive case: 0o = '^'^O) ''^o being a real number in the interval [0, 2]. 
In the attractive case the saddle point is a real number because the resulting 
Hubbard-Stratonovich transformation will not involve the imaginary unity i 
multiplying the linear term in as in the preceding equations. The saddle 
point equation is given by 

p „ J-oo %uJn + n — e — UnQ 
Eq. (|5.8| ) is just a Hartree-Fock aproximation. This means that in the case of 



an infinite range interaction of the form described in Eq.(5^) the Hartree- 
Fock approximation is the exact solution of the problem. Note that in con- 
trast to the case of the HK model, we have a Fermi liquid in the present 
situation. 



6 Conclusion 

Strongly correlated electron systems are very important in today condensed 
matter physics. Unfortunately, even the simplest models developed in this 
field are very complicated to handle. For this reason, many important fea- 
tures cannot be understood exactly. For example, the MIT in the Hubbard 
model can be studied only approximately. The exactly soluble one dimen- 
sional Hubbard model does not exibits a MIT at half-filling. However, the 
MIT in the Hubbard model can be simulated through the study of exactly 
solvable toy models like the one we have discussed in this paper, the HK 
model. We discussed the HK model within the framework of the Green func- 
tion and functional integration formalism. The metal-insulator transition at 
half-filling is easily characterized from the structure of the poles of the Green 
function. The exact Green function is used for evaluate the static suscep- 
tibility at the nesting wave vector at zero temperature. It was found that 
the susceptibility diverges as U approaches the critical value Uc, signahng a 
localization transition. This also happens in the Brinkman-Rice approxima- 
tion in the Hubbard model. However, in the Brinkman-Rice case the static 
susceptibility is evaluated at zero wave vector rather than in the nesting vec- 



12 



tor. Also, in the case treated here we do not have a power law behavior of 
the susceptibility. 

Another interesting result is that the HK model with infinite range hop- 
ping is equivalent to the Hubbard model with infinite range hopping, giving 
the same physical information. For example, the free energy of both models 
coincide. Actually, this is an expected result once the IRH Hubbard model 
separate the k = modes of the k^O modes and the HK model merges the 
exact atomic and band limits of the Hubbard model. Thus, we have shown 
that the HK model shares one more propertie with the Hubbard model. In 
the limit of unrestricted hopping the two models arc equivalent. 

Finally, the HK model admits a further generalization. It can be solved 
exactly even in the case of an arbitrary number of components. 
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